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Abstract 



Compact U(l) lattice gauge theory in (2+1) dimensions is studied on 
anisotropic lattices using Standard Path Integral Monte Carlo techniques. 
We extract the static quark potential and the string tension from 1.0 < At < 
0.333 simulations at 1.0 < /3 < 3.0. Estimating the actual value of the 
renormalization constant, (c = 44), we observe the evidence of scaling in the 
string tension for 1.4142 < /? < 2.5; with the asymptotic behaviour in the 
large-/? limit given by K^/J3 = . Extrapolations are made to the 

extreme anisotropic or "Hamiltonian" limit, and comparisons are made with 
previous estimates obtained by various other methods in the Hamiltonian 
formulation. 
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I. INTRODUCTION 



Classical Monte Carlo simulations of the path integral in Euclidean lattice gauge 
theory have been very successful, and this is currently the preferred method for ab initio 
calculations in quantum chromodynamics (QCD) in the low energy regime. Monte Carlo 
approaches to the Hamiltonian version of QCD propounded by Kogut and Susskind have 
been less successful, however, and lag at least ten years behind the Euclidean calculations. 
Our aim in this paper is to see whether useful results can be obtained for the Hamiltonian 
version by using the standard Euclidean Monte Carlo methods for anisotropic lattices [Q, 
and extrapolating to the Hamiltonian limit in which the time variable becomes continuous, 
i.e. the lattice spacing in the time direction goes to zero. The Hamiltonian version of lattice 
gauge theory is less popular than the Euclidean version, but is still worthy of study. It can 
provide a valuable check of the universality of the Euclidean results f^, and it allows the 
application of many techniques imported from quantum many-body theory and condensed 



matter physics, such as strong coupling expansions [^, the t-expansion [Q, the coupled- 
cluster method and more recently the density matrix renormalization group |^ (DMRG). 
None of these techniques has proved as useful as Monte Carlo in (3+1) dimensions; but in 
lower dimensions they are more competitive. 

A number of Quantum Monte Carlo methods have been applied to Hamiltonian lattice 
gauge theory in the past, with somewhat mixed results. A "Projector Monte Carlo" approach 
TO| , [TT|| using a strong coupling representation for the gauge fields and a Greens Function 



Monte Carlo (GFMC) approach |T3|,|Tj|, which makes use of a weak coupling representation 
of the gauge fields. However both the approaches run into difficulties in that the former 
approach runs into difficulties for non-Abelian models, and the later requires the use of a 
"trial wave function" to guide random walkers in the ensemble towards the preferred regions 



of configuration space |T5[ . This introduces a variational element into the procedure, in that 
the results may exhibit a systematic dependence on the trial wave function |lT6|-pO[]. For this 
reason, we are forced to look yet again for an alternative approach. 



2 



As mentioned above our aim in this paper is to use standard Euclidean path integral 
Monte Carlo (PIMC) techniques for anisotropic lattices, and see whether useful results can 
be obtained in the Hamiltonian limit. Morningstar and Peardon [^] showed some time ago 
that the use of anisotropic lattices can be advantageous in any case, particularly for the 
measurement of glueball masses. We use a number of their techniques in what follows. 

As a first trial of this approach, we treat the U(l) gauge model in (2+l)D, which is one 
of the simplest models with dynamical gauge degrees of freedom, and has also been studied 
extensively by other means (see Section II). Path integral Monte Carlo (PIMC) methods were 
applied to this model a long time ago by Hey and collaborators [^^ , but the techniques 
used at that time were not very sophisticated, and the results were rather qualitative. Very 
little has been done since then using this approach on this particular model. 

The rest of this paper is organised as follows. In section II we discuss the U(l) model 
in (2+1) dimensions in its lattice formulation, and outline some of the work done on it 
previously. The details of the simulations, including the generation of the gauge-field config- 
urations, the construction of the Wilson loop operators, and the extraction of the potential 
and string tension estimates are described in section III. In section IV we present our main 
results for the mean plaquette, static quark potential, string tension. The static quark po- 
tential has not previously been exhibited for this model, as far as we are aware. Finally we 
make an extrapolation to the Hamiltonian limit, and comparisons are made with estimates 
obtained by other means in that limit. (Sec. IV). We find that indeed the PIMC method 
can give better results than other methods, even in the Hamiltonian limit. Our conclusions 
are summarized in Sec. VI. 

II. THE U(l) MODEL 

Consider the isotropic Abelian U(l) lattice gauge theory in three dimensions. The theory 
is defined by the action 
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S = /3Y1 ReTrP^. (1) 



r,fj,,u 

where 



1 - ReTr{\J^{r)\J,{r + fi)\jUr + z>)UUr)|] (2) 



is the plaquette variable given by the product of the hnk variables taken around an elemen- 
tary plaquette. The hnk variable U^(r) is defined by 

U/,(r) = exp[zeaA^(r)] = exp[i6'^(r)] (3) 

where in the compact form of the model, 6^{r){= eaA^{r)) G [0, 27r] represents the gauge 
field on the directed link r —>■ r + fi. The parameter /5 is related to the bare gauge coupling 
by 

/5 = A (4) 
9 

where = ae^, in (2+1) dimensions. 

The lattice U(l) model in (2+1) dimensions has been studied by many authors, and 
possesses some important similarities with QCD (for a more extensive review, see for example 
ref. Ijl2|)- If one takes the "naive" continuum limit at a fixed energy scale, one regains the 
simple continuum theory of non-interacting photons [^]; but if one renormalizes or rescales 



in the standard way so as to maintain the mass gap constant, then one obtains a confining 
theory of free massive bosons. Polyakov showed that a linear potential appears between 
static charges due to instantons in the lattice theory; and Gopfert and Mack ||2^ proved that 
in the continuum limit the theory converses to a scalar free field theory of massive bosons. 
They found that in that limit the mass gap behaves as 



/87r2 / TT^ , . 

aniD = W^exp -v{0) (5) 



^^^^ . 
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while the string tension is bounded by 



9^ vr^ 



where v{0) is the Coulomb potential at zero distance, and has a value in lattice units 

f (0) = 0.2527 (7) 

for the isotropic case. They argue that eq. (H) represents the true asymptotic behaviour of 
the string tension, where the constant c is equal to 8 in classical approximation. The theory 
has a non-vanishing string tension for arbitrary large (3, similiar to behaviour expected for 
the string tension in non-abelian lattice gauge theories in 4-dimensions. 
In the small-/3 limit the string tension behaves as: 

aV = -ln/3/2 + .... (8) 

and the large-/? behaviour is given by eq. (^. 

Similar behaviour will apply in the anisotropic case. Generalizing discussions by Banks 
et al and Ben-Menahem [^], we find (see the appendix A) that the exponential factor 
takes exactly the same form in the anisotropic case, with a Coulomb potential f (0) changing 
to 0.3214 in the Hamiltonian limit. 



For an anisotropic lattice, the gauge action becomes [47,48 



S = PsY.M^) + Y.PtPu{r) (9) 

r,i<j r,i 

where Pij and Pa are the spatial and temporal plaquette variables respectively. In the 
classical limit 

/3s = = —At 

A = ^ = (10) 

e^ttt At 

where At = at/ as is the anisotropy parameter, is the lattice spacing in the space direction, 
and at is the temporal spacing. The above action can be written as 



S = (3 



^^E E ( 1 - cos%(r) ) + ^ ^ ( 1 - cosOuir] 



(11) 



In the limit At 0, the time variable becomes continuous, and we obtain the Hamilto- 
nian limit of the model (modulo a Wick rotation back to Minkowski space) . The Hamiltonian 



version of the model has been studied by many methods: some recent studies include series 
expansions fS^J^, finite-lattice techniques ||30|, the t-expansion pT|j3^ , and coupled-cluster 
techniques [^,|^, as well as Quantum Monte Carlo methods p3| , |36| , |37| , p!2| . Quite accurate 
estimates have been obtained for the string tension and the mass gaps, which can be used 
as comparison for our present results. The finite-size scaling properties of the model can be 
predicted using an effective Lagrangian approach combined with a weak-coupling expansion 
|38| , and the predictions agree very well with finite-lattice data [l^ . 



Some analytic results for the model in the weak-coupling approximation appear in the 
appendix. 



III. METHODS 
A. Path Integral Monte Carlo algorithm 

We perform standard path integral Monte Carlo simulations on a finite lattice of size 
A^^ X Nt, where Ng is the number of lattice sites in the space direction and A^^ in the 
temporal direction, with spacing ratio At = at/ as. By varying Ar it is possible to change 
at, while keeping the spacing in the spatial direction fixed. The simulations were performed 
on lattices with Ns = 16 sites in each of the two spatial directions and Nt = IQ — 64 in the 
temporal direction for a range of couplings (3=1 — 3. 

The ensembles of field configurations were generated by using a Metropolis algorithm. 
Starting from an arbitrary initial configuration of link angles, we successively update link 
angles ^^(r, r) at positions {r,T) which are chosen randomly each time. We propose a 
change to this link angle, which is randomly drawn from a uniform distribution on 
[—A, A], where A is adjusted for each set of parameters to give an acceptable "hit rate" 
around 70-80%. The change is accepted or rejected according to the standard Metropolis 
procedure. 

For high anisotropy (Ar << 1), any change in a time-like plaquette will produce a 
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large change in the action, whereas changes to the space-hke plaquettes will cause a much 
smaller change in the action. This makes the system very "stiff' against variations in 
the time-like plaquettes, and therefore very slow to equilibrate, with long autocorrelation 



times. To alleviate this problem, we used a Fourier update procedure p9| , ^ . Here proposed 
changes are made to space-like hnks which are designed to alter space-like plaquette values 
much more than the time-like plaquette values. At randomly chosen intervals and random 
locations, we propose a non-local change A6'(r, r) = X sin A;(r — tq) on a "ladder" of space- 
like links extending half a wavelength (A = 2TT/k) in the time direction, where both k and 
X are randomly chosen at each update from uniform distributions in [0,7r/2]. We replaced 
approximately 30% of the ordinary Metropolis updates with Fourier updates for anisotropy 
At > 0.444 and 50% for highly anisotropic cases (Ar < 0.444). These moves satisfy the 
requirements of detailed balance and ergodicity for the algorithm. 

A single sweep involves attempting changes to randomly chosen links of the lattice, 
where A^(= SNfN^) is the total number of links on the lattice. The first several thousand 
sweeps are discarded to allow the system to relax to equilibrium. Fig. |I] shows measurements 
of the mean plaquette for /5 = 2.0 and Ar = 1.0, and it can be seen that equilibrium is 
reached after about 50,000 thousand sweeps, with the measurements fluctuating about the 
equilibrum value thereafter. For highly anisotropic cases the system was much slower to 
equilibrate, despite the Fourier acceleration, and in the worst case the equilibration time 
was of order 100,000 thousand sweeps. 

After discarding the initial sweeps, the configurations were stored every 250 sweeps there- 
after for later analysis. Ensembles of about 1,000 configurations were stored to measure the 
static quark potential and glueball masses at each (3 for Ar > 0.4, and 1,400 configurations 
for Ar < 0.333. Measurements made on these stored configurations were grouped into 5 
blocks, and then the mean and standard deviation of the final quantities were estimated 
by averaging over the "block averages", treated as independent measurements. Each block 
average thus comprised 50,000 - 70,000 sweeps. 
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B. Interquark Potential 



The static quark-antiquark potential, V{r) for various spatial separations r is extracted 
from the expectation values of the Wilson loops. The timelike Wilson loops are expected to 
behave as: 

W{r, t) ~ Z(r)exp[ — rV(r)] + (excited state contributions) (12) 

We have averaged only over loops (xq, tq) — > (xq + r, tq) — (xq + r, tq + r) — > (xq, tq + t) — >■ 
(xo, To) which follow either two sides of a rectangle between xq and xq + r or a single-step 
'staircase' route, to estimate W{r,T). To suppress the excited state contributions, a simple 
APE smearing technique |4l|,^^ was used on the space-like variables. In this technique 
an iterative smearing procedure is used to construct Wilson loop (and glueball) operators 
with a very high degree of overlap with the lowest-lying state. In our single-link smoothing 
procedure, we replace every space-like link variable by 



(13) 



where the sum over "s" refers to the "staples" , or 3-link paths bracketing the given link on 
either side in the spatial plane, and P denotes a projection onto the group U(l), achieved 
by renormalizing the magnitude to unity. We used a smearing parameter a = 0.7 and up to 
10 iterations of the smearing process. 

To further reduce the statistical errors, the timelike Wilson loops were constructed from 
" thermally averaged" temporal links That is, the temporal links Ut in each Wilson 



loop were replaced by their thermal averages 

Ut = j dUUe^^{-S[U])/ J dUexp{-S[U]) (14) 

where the integration is done over the one link only, and depends on the neighbouring 
links. For the U(l) model, the result can easily be computed in terms of Bessel functions 
involving the "staples" adjacent to the link in question. This was done for all temporal 
links except those adjacent to the spatial legs of the loop, which are not "independent" p3| . 
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The procedure has a dramatic effect in reducing the statistical noise, by up to an order of 
magnitude 10, which corresponds to a factor 100 in the statistics 

Fig. 1^ shows an example of the exponential decrease of the 4 x 1 Wilson loop W{f, r) with 
r. Estimates of the potential can now be found from the ratios of successive loop values: 



V{r) = i-ln 



W{v,T + At) 



(15) 



W{v,t) 

This ratio is expected to be independent of r for r >> 0. Fig.^ shows the effective potential 
as a function of Euclidean time, r, obtained from Ar = 1.0 at /3 = 2.0 for r = 4 simulation. 
We find that ratio reaches its plateau behavior at very early Euclidean times, which reflects 
a good optimized smearing. We used r = 1 to make our estimates of V^(r). 

IV. SIMULATION RESULTS AT FINITE TEMPERATURE 

Simulations were carried out for lattices of A^^ x A^^^ sites, with A^^ = 16 and A^^ ranging 
from 16 to 48 sites, with periodic boundary conditions. Each run involved 250,000 sweeps 
(350,000 for high anisotropy) of the lattice, with 60,000 sweeps (100,000 high anisotropy) 
discarded to allow for equilibrum, and configurations recorded every 250 sweeps thereafter. 
Coupling values from (3 = 1.0 to 3.0 were explored at anisotropics Ar ranging from 1 to 
1/3. We fixed Ar = 16/A^( in the first pass, so that the lattice size T = Ntat in the time 
direction remains constant, corresponding to a fixed (low) temperature. 



A. Mean Plaquette 

The mean plaquette, or equivalently the 1x1 Wilson loop, was measured for various 
16 X 16 X Nf (with A^^ = 16 up to 48) sizes and for a variety of (3 values, ranging from strong 
coupling to weak coupling regime. Fig. ^ shows the behaviour of the mean spatial plaquette 
< P > for different (3, at fixed Ar = 1 (isotropic case). At small coupling (large (3) the 
plaquette expectation value is expected to approach 1, and conversely it approaches in 
the high coupling (small (3) limit. The plot also shows the 0{f3'^) strong-coupling fSlI] and 
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0(l//9^) weak coupling expansion. The dashed curve represents the weak-couphng series 



and strong-couphng series is represented by the sohd curve. As can be seen that the data 
follow closely the weak coupling expansion down to (3 = 2.0 and strong coupling expansion 
upto (3 = 1.70 and the results agree very well. The variation of < P > with coupling is 
extremely smooth, with no sign of phase transition, as we should expect. The cross-over 
from strong to weak coupling seems, from the Monte Carlo data, to take place quite raidly in 
the region /? ^ 1.8 — 2.0. Horsley and Wolff |^ investigates the effect of finite size lattice in 



their weak-coupling expansion calculations. They find that such effects enter at order l//3^ 
as a correction of order l/l^, where I is the lattice size and D is the number of dimensions. 
Thus for U(l) in (2-|-l) and / > 10, finite size effects should be essentially negligible. 
Fig. ^ shows a plot of our estimates of < P > as a function of anisotropy Ar^ for the case 
P = We would like to make contact with previous Hamiltonian studies by showing that 
the mean plaquette value approaches to previously known values in the Hamiltonian limit 
Ar 0. The extrapolation was performed using a simple cubic fit in powers of Ar^. In 
this limit our results agree very well with the Hamiltonian estimate obtained by Hamer et al 
[r?! . The estimates of the mean plaquette < P > in the Hamiltonian limit are graphed as a 
function of /5 in Fig.^ and compared with the series estimates at strong and weak coupling 
28| , and Greens function Monte Carlo results |T^. At large /3, the data is in good agreement 



with predictions of weak-coupling perturbation theory. At small /3, the data follow closely 
the strong-coupling and Greens function Monte Carlo estimates, except at very small /3. 
We attribute the discrepnacy to the finite-temperature effects and a better agreement is 
expected at zero-temperature [52]. 

B. Static quark potential and string tension 

The effective potential is obtained from the ratio of successive Wilson loop values: 



/Vr(r,r + Ar) 
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Fig. 1^ shows a graph of the static quark potential V(R) as a function of radius R at /3 = 2.0 
and At = 1.0. To extract the string tension, the curve is well fitted by a form 

V{R) = a + b\nR + aR, (17) 

including a logarithmic Coulomb term as expected for classical QED in (2+1) dimensions 
which dominates the behaviour at small distances, and a linear term as predicted by Polyakov 
26| and Gopfert and Mack [^] dominating the behaviour at large distances. The linear 



behaviour at large distances is very clear, but the data do not extend to very small distances, 
so there is no real test of the presumed logarithmic behaviour in this regime. 

Figs. 1^ shows the behaviour of the fitted value of the string tension K{= a'^a) as a 
function of P for the isotropic case ( At = 1), plotted in comparison to results predicted for 
string tension using the large-/? approximation (eq. |^). We also include the strong-coupling 
values and show the leading order strong-coupling prediction (eq. 0). As can be seen that 
our estimates of the string tension, in the weak-coupling region are much larger than the 
Villain prediction, though slopes are more or less well determined. However one does not 
expect agreement of precise magnitude, since our estimates are obtained using Wilson ac- 
tion, and not the Villian approximation. But one does expect that the slope of the string 
tension in both cases should ultimately be the same as one approaches the continuum limit. 
Mack and Gopfert argued that eq. (^ represents the true asymptotic behaviour of the 



string tension, apart from the value of the renormalization constant c. Such a behaviour 
was obtained by a semi-classical treatment of the theory with the effective action without 
the correction terms. Thus it is possible that eq. (^, including the normalization constant, 
becomes exact as 1//3 — 0, if one estimates the actual value of c. We estimated the actual 
value of c by graphing the e^'^^"^^ x AAA3K against 1/(3 and extrapolating to 1//3 — > 0. Using 
estimated value of the renormalization constant, c = 44, we plotted the asymptotic form of 
the string tension in weak-coupling limit in comparison to our estimates. The result agrees 
very well with our estimates. 
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Fig. 1^ shows the behaviour of the string tension (K) as a function of the anisotropy Ar^, 
for fixed couphng f3 = V2. An extrapolation to the Hamihonian hmit At ^ is performed 
by a simple cubic fit. Again the extrapolation agrees rather well with earlier Hamiltonian 
estimates [|12|. Note that this quantity depends rather strongly on At: there is a factor 
of three difference between the values at At = and At = 1. Extrapolating to At 0, 
estimates of the string tension in the Hamiltonian limit, for various (3 values are obtained. 
The plot of these estimates, in comparison to the large-/? prediction of eq. (|^) is shown in Fig. 
p!0| . As can be seen, the weak-coupling results, with estimated value for constant coefficient, 
agree very well with our estimates. The comparison with the theoretical prediction using 
Villian form is not bad either, especially in the large-/3 region. In this region, upto (3 ~ 1.8, 
the results agree well within errors and slopes are very well defined. 



Fig. [Tl| shows the plot of Hamiltonian estimates as a function of (3, plotted in comparison 
to the previous estimates obtained from strong-coupling [8], Greens function Monte Carlo 



| r?| and finite-size scaling The agreement with the series expansion and GFMC results 
is once again excellent. The estimates of finite-size scaling describe the data at /? = 2.0. 



V. SUMMARY AND CONCLUSION 

In this paper, we used numerical simulation techniques on anisotropic lattices for U(l) 
model in (2-1-1) dimensions to obtain a clear picture of the confinement and improve our 
knowledge of the string tension. Our first aim was to implement Path Integral Monte Carlo 
method to the U(l) model, which is one of the simplest models and has been studied by 
various other methods. It seems to work extremely well for the mean plaquette in the 
strong-coupling and weak-coupling regimes. The simulation results were extrapolated to 
At 0, and Hamiltonian estimates for mean plaquette were obtained. A comparison was 
made with the previous estimates obtained in this limit and agreement was excellant. The 
second aim was the get a better picture of the confinement and the numerical study of the 
asymptotic behaviour of the string tension in the large-/? limit. It was shown that at large 
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separations the U(l) theory in (2+1) dimensions shows a hnear behaviour as predicted by 
Polyakov. However, the presumed logarithmic term could not be tested as our data do not 
extend to very small separations. The estimates of string tension have confirmed very well 
the weak-coupling predictions, provided that one takes the estimated value of the constant 
coefficient c rather than its value in the classical approximation using effective action without 
correction terms. For /3 < 2.0, the data is in good agreement with strong-coupling 
predictions. Once again our Hamiltonian estimates for the string tension fit very well with 
the series and GFMC results. 

This seems to suggest that PIMC does work better than other methods, as stated by 
Ceperley. The first clear picture of the static quark potential has been obtained, showing 
clear evidence of the linear confining potential at large distance predicted by Polyakov. Our 
results from the PIMC do appear to be in agreement with results of the series expansions. 
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APPENDIX A: WEAK-COUPLING PREDICTION USING VILLIAN 

APPROXIMATION 



Here we extend the discussion of the weak-couphng behaviour of the mass gap to the 
anisotropic case, following the discussions of Banks, Myerson and Kogut |^ and Ben- 
Menahem |^^. This can be probed by means of the expectation value of the Wilson loop, 
< W{C) >, which for pure gauge theory is given by 

Z{J) 



< W{C) >-- 



where 



Z{J)=l[{ fd0,{r)]exp 



m 



r,iJ. 



(Al) 



(A2) 



and J^{r) is a vector field associated with a rectangular contour C of width R and length 
T. The extra piece of action associated with the external charge is 

S.J = I d'rJJr)eJr) 



For a point charge, Sj is given by the path integral of the gauge field along the world line 
of the charge. Therefore we choose 



Mr) 



1 if r — i> r + /i is on C 
— 1 if r + /t ^ is on C 
otherwise 



(A3) 



and this current is conserved: A^J^(r) = 0. We will assume a spacelike Wilson loop for this 
discussion. 

It is convenient to work in the temporal gauge 



9t(f^ = 



(A4) 



Then one can separate the angular variable 6i{r) into longitudinal and transverse parts 

1 
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(A5) 



whence it follows that = 6{r). In the temporal gauge eq. ( [X^ ) can be written as 



Z(J) = n/ Wr)}exp 



-pAr ^ ( 1 - cos 



2Ar V"" 'A,A, 



(A6) 



with = + A2. A and A are defined are the lattice difference operators. To calculate 
the path integral we replace the potential term by its Villain form 



-/3(l-cose(r)) 



(A7) 



l=—oo 



where is the modified Bessel function and exp(i/^) is the character of the plaquette 

variable in the Ith irreducible representation of the compact U(l) group. 
For large (3 the modified Bessel function can be approximated by 

1 



(A8) 



The Gaussian integration over 9{r) reduces the path integral to 



z{j) oc n E 

/{r)=— 00 



2(3 At ^ ^ \ At At J ^ ' P r AiAt ^ 



(A9) 



The convergence of the sum over / can be improved by making use of the Poisson summation 
formula, 



00 00 „ 

Y: h{l) = E y d<ph{cP)e 



2iTTmr 



(AlO) 



m=— 00 



l=—oo 

where h is an arbitrary function. 

Substituting /(r) = AtL{r) and using eq. (|A8|) , Eq. (|A9D becomes 

1 



/oo ^ 

' m(r)=— 00 



exp 



2i3At 



J2<l>{r){AT^A,A, + AtAt](t>(r) 



xexp 



(All) 



Finally doing the Gaussian integral over yields; 

r 

1 



z{j) (X n 

m{r)=—oo 



1 



+27riAr^Zm( 



1 



m{r) 



Ar^AAr + ^AtJ^'"'^'^' 
1 



xexp 



-—yj- 



2/5 'Ar^A.A. + AiAi 
The partition function, Z( J = 0), is given by 



Mr) 



(A12) 



where 



m(r)=— oo 



^(r - r') 



-27r2/? ^m(r)G'(r - r')m(r') 



(A13) 



Ar^AiAi + Ai At 

is a three-dimensional lattice massless propagator. 

The propagator can be expressed as a momentum integral: 

'■^ dkQ r dh r /^^^Qik.r+kot) 



(A14) 



J-T7 ZTT J-77 2/K J-TX 



4 sin^ A;o/2 + Ar^ ^4-2 cos k^-l cos /ca) 



(A15) 



For the isotropic lattice (At = 1), G{k = 0) = 0.2527; while in the hmit At oo, G{k = 0) 
takes the value 0.3214. 

The resulting partition function may be written as 



Z(^J) ^spin—wave^ext^mcmopoles 



(A16) 



where 



^ spin— wave 



j[dx] 
Zext = exp 



exp 



Y.Ji{r)G{r-r')Ji{r') 

' r.r'i 



Zmon^ole^Y. E ^Xp 

m(r)=— oo 



-27r^/3 ^ m{r)G{r - r')m{r') + 2m ^ G{r - r')r]{r)m{r) 



(A17) 
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The first factor in eq. ( |A15| ), tlie "spin- wave" piece, contributes tlie ordinary Coulomb part 
to tlie potential between the external charges. The rest of the eq. ( [A 151 ) is the partition 
sum of the Coulomb gas of magnetic monopoles interacting with a stationary current loo 
via G{r — r ). This three-dimensional Coulomb gas is always in a plasma phase. The only 
thing that changes with /3 is the density of the monopoles which decays exponentially to 
zero as (3 increases. The monopole density generated by the external loop and coupling to 
the lattice monopoles m(r) is a dipole sheet along the surface spanning the contour C. Thus 
eq. ( |A15|) represents a dipole sheet immersed in a monopole gas, which interacts directly 
with the dipole sheet. The gas becomes polarized, and the monopoles of opposite polarity 
accumulate on the two sides of the sheet, screening both the magnetic field away from the 
sheet and the monopole-monopole interaction, thus minimizing the energy of the system. 
This gives the correlation function a finite range, whose inverse is interpreted as a mass of 
the scalar field. 

In the naive continuum limit, the Coulomb self-energy of the monopoles blows up and their 
density goes to zero exponentially. This causes the coefficient of the linear force law to go 
to zero and one regains the usual continuum theory, i.e., free electromagnetism in three 
dimensions; but if one renormalizes or rescales in the standard way so as to maintain the 
mass gap constant, then one obtains a confining theory of free massive bosons. The detailed 



analysis of Polyakov [|26[ shows that the sheet cannot be screened completely and there is 
a finite energy left that binds the quarks through a potential which grows linearly with 
separation. Polyakov showed that for arbitrarily large finite (3, theory has a mass gap and 
there appears a linear potential appears between static charges due to instantons in the 
lattice theory. 

Splitting the exponent in eq. ( |A15| ) into pieces with r = r' and r 7^ r', the monopole 
partition sum can be written as 



■'monopole \ 



(J)=n/#Wexp - ^ A,0(r) - Y: ^^'PGo 
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r)m[r 



r m r 



(A18) 



For large /?, one can neglect all the monopoles with charge other than 0, ±1 in the sum, 
since the monopoles with a higher charge gets extra powers of e~^. Q Applying the large-/? 
approximation to the monopole partition function, the above equation can be written, in 
the naive continuum limit, as 

-,2 



(J) ^ j D(Pexp 1 1^(9^0)2 -M2 cos (^0 + 27rr/)|c?3x 



^monopole \ 

where the lattice photon mass M is given by 



(A19) 



2 _ StT^ ^-2n^Go/g^a_ 



(A20) 



At this point out analysis parallels the work of original work of Banks et al 
The the expectation value of the Wilson loop becomes 



< W{C) > ^ exp 



1"^ r.r'i 



Y[ / (i0(r)exp 



2g-2^2/3Go ^ ( 1 - cos(0(r) + 27Tr]{r)) 



(A21) 



APPENDIX B: LARGE-/? PREDICTION IN HAMILTONIAN FORMULATION 



We now do the Hamiltonian study of U(l) model in (2+1) dimensions using weak coupling 
perturbation theory. The Hamiltonian formulation of compact U(l) gauge theory has been 



studied by Banks et al. |^5[ and Drell et al. pj] following the Kogut and Susskind approach 



Consider the anisotropic Wilson gauge action for U(l) model in (2+1) dimensions, 

1 



s = p 



At 



(Bl) 



^ This argument becomes problematical when Ar ^0 [46|. In this hmit the monopoles of higher 
charges contribute to the action and hence are no longer suppressed. 
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where we made assumed the lattice symmetric in spatial directions with = 1 and at = At. 
The plaquette variable is given by 



Pij = 1 -cos^i; 



and f3 = l/g\ 

In the weak coupling approximation 



cos^5^1-i^f 



The plaquette angles 6^ are related to the link angles Ai by 



9^{n', k) = eijk Ai{u,i) - Ai{n + j,i) + Ai{n + i,j)Ai{n,j) 



(B2) 



(B3) 



where (n', k) denote a plaquette centered at n and oriented in the k while {n ,i) describes 
a link that starts at site h in the direction i, with 



n = n + 



.i + j- 



The Fourier transformations of the plaquette angles and the link angle give 

1 



0-M 



where N = N^Nf is the total number of sites; then 



Alii,) = 



Therefore eq. ( |B3[ ) can be written as 



where 



(B4) 



(B5) 



(B6) 
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y)^2iY, e,.,e-'^'^'^ sm(^,/2) 



(B7) 



Choosing the temporal gauge Aq — 0, the time-hke plaquette contributions, in the weak 
couphng approximation, become: 



^4(i)A^(i)x4sin^(W2) 



(B8) 



while the space-like plaquette contribution becomes: 



" k 

= E4^^^^ 



(B9) 



where and Rj, are the vectors and matrices, (2 x 2) in (i,2) square, respectively, and 



i?- = gIG-^. Therefore 



( 



shi^{ky/2) _ei(fe./2-fc,/2) ^^12 sin ky/2 

.Qi{ky/2-k^/2) g-^ j^^i2 gin ky/2 sin^{k:,/2) 

The matrix R^^ is hermitian, and has eigenvalues: 

A^(l)=4X:sin^(W2) 



(BIO) 



(Bll) 



We can define a unitary transformation 



= ^k^k 



(B12) 



which diagonalizes the matrix i?^: 



_eifex/2gin(^.^/2) e^'=^/2)gin^^/2 \ 



sgrisin^ kx/2 + sin^ /cj^/2 
The components of ujy. are 



V 



.Qik^m gij^ ^^^2 6*^=^/2 sin2(A;^/2) 



(B13) 
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i 



Cl 



e"*'^-/') sin KI2A^(2) - e'^-'^ sin(A;^/2)A^(2) 



where Ci = ysin^ A;^/2 + sin^ 

The component of wj, corresponding to the zero eigenvalue is 



Cl 



This vanishes at all orders and is an "ignorable" coordinate. 
Thus the action eq. ( |Blj ) can be written as 



5 = fE[Ar(A^^(lK(l) 



+ 



A 



r 



4sin2(fco/2)^u;|(/i)a.^(/i) 



+ const. 



where 



sin^ A;i/2 + sin^ ^2/2 



Therefore the partition function can be written as 

Z = \{{j duk]e-^^^^ 

k 

where the action can be written in the following form: 



2 ^ " ^ 



+4(2V^,(2)(-^sin^(V2; 



4 

Now ujj, is the ignorable coordinate and will be ignored from here on. 
The partition function becomes 



k 



exp 



27r 



1 



/5 



^^^k + i?sin2(V2) 



1/2 



f In(^) - ^El^((^^Ak + ^sin2(A:o/2)^ 
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This lead to the following free energy density: 



27r 



d{k) 



(B21) 



where 



d{k) = (ArA^ + — sin2(A;o/2) 



(B22) 



1. Mean space- like plaquette 

At weak coupling the cosine term in the space-like part may be expanded in the powers 
of ^p. Therefore 

1 9 

< cos % > = 1 - - < dtj > 

k 

where 

< 4 i Vfe 1 > = ^ . t ^ B24 

' ' /du;^e-t''('=H-^ 

Using the integrals 

due-"^ =J- 
-<x> V 2 



we get 



< cotu;-, >= ^ (B26) 



pd{k 



Therefore 



< cos >= 1 - y (B27) 

Now consider the "cylindrical limit", Nt oo, At — > 0: In this limit the sum term in the 
right hand side of the above equation becomes; 
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A- ^ /■2t 



d{k) Y 277 7o "Ar + ^sm^(V2) 



^0 (ArA(^) + £) 



(B28) 



In the bulk limit, N oo, this becomes: 



^ 2VA(A:) 



y (ikj-dky J 4(^8111^ kx/2 + sin^ ky/2 



I 

2N'^C["> (B29) 



2 



where 



2Ar2 
* k 



0.958 (B30) 



is a standard lattice sum. 
Therefore 

< cos^i,- > = 1 - -^Ci^^ 
^ 2(3 ^ 

_ 1 3 ^(2) 

= 1 - 0.479^^ (B31) 

This result agrees with the estimates in the Hamiltonian limit Now taking the bulk limit 
in isotropic case, we get 

^^^J^iridkdk d- (sin^ fc./2 + sin^ fc,/2) 

^ dQz) (27r)3iio i ^ "(sin2A;^/2 + sin2A;y/2 + sin2A;o/2) ^ ' 



By symmetry, we can easily find that 
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Therefore 



KcosOij >^1-^ (B34) 



2. Mean time-like plaquette 



In weak coupling approximation the time-hke plaquette can be written as 

1 2 1 2 

- ^ < COS^iO > = 1 - T E < ^'o > 
^ i=l ^ i=l 

= l-i<^E4sin^W24(^M^(^)> 

it 



1 - 4^ < E^sin^ W2(c^i(iK(i) +a;i(2K(2) > (B35) 



Now 



< 4(1)^^(1) + 4(2)c^^(2 >= —rr + 



At 



I3d{k) 4/5sin^A;o/2 



(B36) 



^ E < cos^io >= 1 - -^E^sin' W2 



2 

47V 



»=i " it 

The term in the square brackets can be written as 

N r r^^ 



1 At 



E4sin2A;o/2 

k 



1 At 

+ 



_d{k) 4sin'A;o/2 



(27r)-J JO 



4sin'A;o/2. 

dk^dkydkoA sin^ /cq 
1 



(B37) 



sin ko/2 



sin^ A;o/2 + sin^ ko/2 + AT^(sin2 k^/2 + sin^ ky/2)_ 
(At = 0) 

N{l + l/3) (Ar = l) 



(B38) 



Thus the mean time-like plaquette gives; 



1 2 

-E < cos 61^0 >= < 
^ i=i 



(At = 0) 
1-^ (Ar = l) 



(B39) 



as is expected. 
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3. Expectation value of Wilson loop 



The expectation value of the Wilson loop is given by 



(B40) 



where the Wilson loop operator is given by 

1 



W, 



c 



n e^' + h.c 

i&c 



(B41) 



Taking the separation R in the x direction, for simplicity, the Fourier transformation of 
above equation) is 



Wc = exp 



+ h.c 



(B42) 



k x,^iec 

We now consider the contribution from class of "zero modes" with k = {kQ,k = 0). These 
modes are likely to be responsible for the dominant finite-size corrections 0, i.e.. 



Wo{R,t) = exp 



EA-o,o(^)^h 



ko 



Jkoxoi „ikoxo2 



(B43) 



Recalling that Aq = 0, the above equation becomes 

2 



Wo{R,t) = exp 



exp 



ko 



5]wfc(i)i?sin koT/2 



ko 



(B44) 



where T = xq2 — xqi and we have ignored the e 2 >. Therefore 



<WoiR,T) >-- 



^ny^-^-o(i)exp 



uJkoRsm{koT/2) - ^^Asm^{ko/2)ujl^uJko 



(B45) 



Let 



2At 



4sin2(fco/2) = c{ko); 



then the numerator of above equation can be written as 
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duexp 



-c{ko (^uj — i 



Rsm{koT/2)' 







)] 


exp 



sm\koT/2) 



Nciko) 



Therefore the expectation value of the Wilson loop is given by 

-n- ri?2sin2(A;oT/2)' 
<Wo{R,T)> =nexp ^ ° ^ ^ 

fco 

exp 



N2l3sm\ko)_ 
R^At Nt /•2^ sin2(A;oT/2 

u/Cq — ; — 2 



2N^Nt(3 27rJo " sm'{ko/2 
Now consider the cylindrical limit, Nt oo, At — > 0. In this limit 

271 Jo " sin2(A;o/2) 

Therefore 



< Wo{R,T) > = exp 



2 'N, 
exp I — (T X area 



TAt 



Where 



^eff 



R_^ 
2 ^N?' 



2 ^L2' 



(B46) 



(B47) 



(B48) 



(B49) 



(B50) 



For R=L, i.e., for pair of Polyakov loops mapping around the lattice in the spatial direction, 
we get 

-,2 



2L 



(B51) 



This corresponds to the value calculated by Hamer and Zheng 0, a = 1/L, for a string 
wrapping around the entire lattice, if we multiply by an extra factor g'^/2 coming from 
renormalizing Hamiltonian. As is obvious that for smaller loops, the effective string tension 
is less. 
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FIG. 1. Plot of plaquette average against the number of configurations. The loop average is 
seen to reach its equilibrum value, after discarding the initial number of sweeps (50000 for 16^ 
lattice) . 



30 




FIG. 2. Plot of a 4 X 1 Wilson loop against r from At = 1.0 simulations for /? = 2.0 
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FIG. 3. Ratio of the Wilson loops as a function of r for fixed x and y from /? = 2.0 and 
At = 1.0 simulations. The ratio is seen reaching its plateau at very early times and is very close 
to its asymptotic large-t plateau. 
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P 

FIG. 4. Plot of average plaquette as a function of (5. The o symbols label our Monte Carlo 
estimates. The solid curve represents the 0{(3^) strong-coupling expansion [49] and the dashed 
curve represents the 0(l//3^) weak-coupling expansion [50]. 
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FIG. 5. The mean plaquette as a function of Ar^ at /3 = 1.4142. The sohd curve is cubic fit, in 
powers of Ar^, to the data extrapolated to the Hamiltonian limit. A symbol represents the series 
estimates [27] in that limit. 
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FIG. 6. Mean plaquette estimates against /3 at At = 0. The earlier estimates from 
strong-coupling expansion and Greens function Monte Carlo are represented by O and □ sym- 
bols respectively. Our Monte Carlo estimates and weak-coupling expansion results, (Eq.BSl), are 
represented by Q ^^'^ ^ respectively. 
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V(R) 




R 

FIG. 7. The static-quark potential, V{R) as a function of the separation R. This plot includes 
the measurements from /? = 2.0 for At =1.0 with 10 sweeps of smearing at smearing parameter 
a = 0.7. The analysis errors are smaller than the plot symbols. The dashed line is the fit to the 
form V{R) = a + bR + dog{R). 
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FIG. 8. String tension as a function of j3. The solid curve represents the weak-couphng result 
using Villian approximation, with our estimated value for renormalization constant, c(=44) and 
the dotted curve coinciding with the solid curve is the linear fit to our Monte Carlo estimates. 
The dashed curve represents the scaling behaviour predicted by Mack and Gopfert (eqj^ and the 
dash-dot curve is the strong coupling result of eq. (jsl). 
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FIG. 9. String tension against Ar^ at (3 = 1.4142. Extrapolation to Hamiltonian limit is 
performed by a cubic fit shown by the dashed line and the earlier Hamiltonian estimate [27] is 
shown by Q- 
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FIG. 10. plot showing the estimates of string tension in the Hamiltonian limit as a function 
of (3. The solid line is the weak-coupling prediction with estimated value of the renormalization 
constant c, and dashed line is the linear fit to our Monte Carlo data. The dash-dot curve represents 
the large-/? prediction of Mack and Gopfert using Villian approximation. 
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FIG. 11. Hamiltonian estimates of string tension as a function of p. Our estimates of string 
tension are labeled by Q symbol. Earlier results from series expansion [8] , Greens function Monte 
Carlo [17], and finite size scaling [37] are labeled by A, □ and ★ respectively. As is seen that our 
estimates make a very good contact with earlier Hamiltonian estimates obtained by many other 
techniques 
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